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The transmitted wave that results from a collision of a wave packet which is initially to the left 
of a potential barrier depends in general on the amplitudes of negative momenta of the initial state. 
The exact form of this dependence is shown and the importance of this classically forbidden effect 
is illustrated with numerical examples. Special care is taken to account properly for bound states. 
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I. INTRODUCTION 



Suppose that a classical ensemble of independent particles in one dimension is initially confined (at t = 0) in the 
spatial interval a < x < b < 0, and allowed to move freely after t = 0. Only particles with positive momenta may arrive 
at positive positions for t > 0. In contrast, the quantum wave function involves negative-momentum contributions as 
well, 



ii(x,t) = h- 1/2 / d P e lpx/h ^{ P )e- lE ? t/h , 



where E p — p 2 /2m, and 



^(p) = h- 1/2 / dxe- ipx/R ip{x,0) 



(1) 



(2) 



is the momentum representation of the initial state. The effect of negative momenta for x, t > is however a transient 
one; the total final probability to find the particle at x > is given only by positive momentum components, 

/>oo />oo 

P T (oo) ee lim / dx \ip(x,t)\ 2 = / dp\ijj(p)\ 2 , 



(3) 

"^JO Jo 

see e.g. M or 0, since there are no bound states. This negative-momentum effect is also present in collisions, where 
it is combined with other classically forbidden effects. Consider the family of cut-off potentials of the form 



r if x < c , 

V{x) = { U(x), if c < x < d, 



(4) 



IVo, 

where Vq > and V(x) are real, as depicted in Fig. |l|. 
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FIG. 1. Scattering process in a "cut-off" potential: the support of the potential is compact; an incident wave with momentum 
p produces a transmitted and a reflected wave. 



1 



Let us suppose that the maximum value of the potential is Vm ■ For a classical ensemble confined between a and b 
(such that a<b<c<d<0) only particles with initial momentum above the "barrier momentum" pm = (StoVm) 1 ^ 2 
may pass to the right of the potential region. In the quantum case however, there are contributions from all the 
components of the initial wave packet: (a) p > pm = (ZjtiVjw) 1 ' 2 (above-the-barrier transmission); (b) po < p < Pm 
("asymptotic" tunneling, these momenta contribute to the transmission probability Pt(oo)); (c) < p < po ("transient 
tunnelling", these momenta do not contribute to Pt(oo)); (d) — oo < p < (transient negative-momentum effect); (e) 
Pj = ijj, 7j > (tunnel effect associated with bound states, which do contribute to Pt(oo)). 

While the various tunnel effects have been well discussed in the literature, even though not as clearly distinguished 
as in the above classification, the negative momentum effect, actually the only one that survives for free motion, has 
been frequently overlooked, a clear exception being . This paper is complementary to the more usual treatments 
of one dimensional scattering |^-^], and lies in the wake of who did consider in a concise manner all these effects, 
with special emphasis on the inclusion of bound states. Our S-matrix treatment is more explicit than that of ||,||, 
and shows how the resolution of the identity in scattering eigenstates leads to compact expressions for the transmitted 
wave packet. We also illustrate the negative momentum effect with some numerical examples and work out in detail 
both the case when the state is confined initially to the lower level and when the initial confining is to the upper level 
(for which case the frequently overlooked contributions are those of positive momenta, including evanescent waves) . 
The contribution of bound states is indicated explicitly, both formally and with some numerical examples. 



II. "STATIONARY" EIGENSTATES OF THE HAMILTONIAN 



The total Hamiltonian H = Hq + V may have a discrete set of bound states {\Ej)}, with energies Ej < and (real) 
wavefunctions <j>j(x), and a continuum of "stationary scattering eigenstates" with E p > 0. Only the former belong to 
the Hilbert space of square integrable functions. The latter however form a convenient basis normalized according to 
Dirac's delta. For energies above Vb, the energy spectrum is doubly degenerate as corresponds physically to incidence 
from one side or the other. Below Vg there is only one linearly independent solution. The resolution of the identity 
may be written in different ways, in particular as 

/—pa r°° r±Po 

dp\p ± )( P ± \+ dp\p ± )(p ± \± dplp^ip^, (5) 
-oo Jp JO 

where the states |p ± ) have energy E p = p 2 /(2m). Continuum and bound states are orthogonal. Moreover, (p ± \p' ± ) = 
S(p — p r ), and (Ej\E'j) = 5jj>. The first two integrals in Eq. (Js[) reflect the double degeneracy signalled above, whereas 
the last integral corresponds to the non-degenerate part of the continuous spectrum. The reason for the sign choice in 
front of the last integral is due to the unified notation we shall now introduce for the generalized eigenstates associated 
with the continuous part of the spectrum. 

The states \p + ) with p > are characterized by an incident plane wave of momentum p from the left. The 
asymptotic behaviour of their wave-functions ip p + (x) is 

, M _ _j_ [ exp(ipx/H) + R l (p) exp(~ipx/h), x < c , fi . 

Vp + w - h i /2 x \ r i(p) exp ( iqx /n), x >d, w 

where q is the momentum with respect to the upper (right) level of the potential, q — (p 2 — 2mVb) 1 / 2 . For the 
evanescent regime, i.e. when p 2 < 2mVb, the plane waves on the right become decaying exponentials so the positive 
imaginary square root is taken. R l (p) and T l (p) are reflection and transmission amplitudes for left incidence. They 
are obtained by solving the stationary Schrodinger equation subject to the boundary conditions specified in (^). 

The states \p + ) for p < —po are defined by having an incident plane wave from the right, oscillating with spatial 
frequency |<?|/?i, where now the negative square root is taken, q — — (\p 2 — 2mVo\) 1 ^ 2 ■ The asymptotic behaviour of 
the corresponding wave-functions is 

,/, +( r ) - J— ( P\ 1/2 x fT r (-p)exp(ipx/h), x<c , . 

WpK) ~h l l 2 \q) \cxv{iqx/h) + R T {-p)cxv{-iqx/K), x > d, [> 

The factor (p/q) 1 ^ 2 is necessary for the proper delta normalization. R r (—p) and T r (—p) are reflection and transmission 
amplitudes for right incidence. Note that the arguments of transmission or reflection amplitudes are always positive 
for states \p + ), independently of the sign of p. 

In both cases, @ and (^), q may be defined as the square root of p 2 — 2mVo with a branch cut that joins the branch 
points p = ±po going slightly below Im(p) = 0. In this way the sign of q is the same as the sign of p for p 2 > p%. 
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The states |p~) are obtained simply by considering p < in expression (^|) and p > po in expression (^). In these 
states there appears an outgoing plane wave instead of an incident plane wave. Thus, the right hand side of Eq. @, 
with p < 0, defines the corresponding states \p~), as those with a left-outgoing plane wave of absolute momentum |p|, 
and can be read as (x), for p < 0. For p < — po, q < 0, whereas for — po < p < 0, q becomes as before a positive 
imaginary number. States \p~) with a right-outgoing plane wave are defined by the right hand side of Eq. (0), with 
p > po, and q > 0. Note that the arguments of the amplitudes R r,l (p) and T r ' l (p) are negative for \p~) states. Since 
the formal boundary conditions that define the states Ip^) are in fact equal, the negative- argument amplitudes will 
be given by the same formal expressions valid for their positive-argument counterparts. For the same reason, we shall 
refer to i? r ' Z (p) and T r ' 1 (p) as "reflection" and "transmission" amplitudes independently of the sign of p, even though, 
on physical grounds, this terminology and notation would only be appropriate for the p > case (only in that case 
does T r ' 1 multiply a transmitted wave and R r ' 1 a reflected wave). 

Reflection and transmission amplitudes are not independent. The unitarity of the S matrix imposes certain relations 
among them. The S matrix elements arc defined as the coefficients multiplying the outgoing plane waves when the 
incident plane wave is normalized to unit flux. When the two channels of the one dimensional scattering are open 
(this happens for p > po), they are given by 

/ . 1 /o 

R l (p) 

S(p) = | ^' " , xl/2 i (M 




(?) ^(P) 



The unitarity of the S matrix, SS^ = 1, implies that 



P\T r (p)\ 2 + \R r (p)\ 2 = l, (9) 

q 

\R l (p)\ 2 + ^\T l (p)\ = l, (10) 

V 

^T r (p)R l ( P y + R r ( P )T l (py = o. (ii) 

For < p < po only one channel is open, and the S matrix reduces to a number, R (p). Unitarity implies in this case 

R l (p)R l (p)* = 1, 0<p<p o . (12) 

All these equations, from (0) to (|l2|), are also valid for negative momenta and relate the amplitudes associated with 
|p~) states. They can also be derived by comparing various Wronskians of the stationary scattering states at different 
regions. 

For the evanescent case, < p < po, one important relation follows by multiplying |p + ) by R l {p)* and using jl^). 
This gives the state | — p~ ) . Equating the coefficients for x > d, 

T l {-p) = T\p)R l {p)\ 0<p<p o . (13) 

Similarly, by taking the complex conjugate of the boundary conditions (^|) and (0) , it is found that ip— p - (x) — ip p + (x)*, 
and comparing the coefficients that multiply the exponentials, 

T r > l (-p)=T r > l (p)*, (14) 
R r \-p)=R r \p)* . (15) 

Another important relation between T l (p) and T r (p) follows by equating the Wronskians of ip p ±(x) and tp_ p ±(x) at 
x < c and x > d, 

T r (p)p = T l (p)q. (16) 

As an illustration of the above, the reflection and transmission amplitudes for the simple step potential (c = d = 0, 
U(x) = 0) are given, for all p, by 

T'(P) = ^; R\p) = P -^: 
q+p q+p 

T r (p) = ^; Rr {p) = l-P. (17) 
p+q p+q 

The reader may easily check relations (j^[f6|) for amplitudes that correspond to the step potential, Eq. (|l7|), as a test 
of their validity. 
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III. A COMPACT EXPRESSION FOR THE TRANSMITTED WAVE FUNCTION 



In this section we shall find, using the relations of the previous one, an expression for ijj{x, t), with t > and x > d, 
assuming that the initial wave function is restricted to a < x < b < c. First we insert the resolution of the identity in 
terms of bound states and scattering states \p + ), 

^x,t) = Y,^( x )( E ^(°» e ~ lEjt/h 

dp4, p+ (x){p+\^(0))e-' lE ^ h 

J 

+ / dp^ p+ (x){p + \^(0))e- lE ^/ h , VxandVi. (18) 
Jo 

Because of the initial restriction of the wave function, the matrix element (p + \ip{0)) may be evaluated with the aid 
ofEqs. (§) and (0), 



(h) V2 ^(p)T r (~p)\ p<- Po . 
$(p) + $(-p)R l (p)*, P >0. 



(p+|V(* = 0)) 



(19) 



Now we use Eqs. (|-@) and (|T|) in @, for x > d, 

V(x,t)-^^(x)(£; J #(o)} e -^ t ^ 



-Pa 



+ h- 1 ' 2 / dp ^ \e lqx ' h + R r (-p)e~ lqx ' h ] T r {-p)*iAp)e~ tE ^ t/h 
q L ' J 



+ h-^ 2 d P T\p)e lqx/h i>(p)+R l {p)*ijj(-p) 
Jo L 

and reorganize the terms as follows, 

i>(x,t) =Y,M*){Ejm = 0))e- lE > t/R 



dp -T \-p)*^{p)e lqx,h e- tE ^ t/n 

-oo q 



e- iB "^ n , x>d, 



(20) 



-Pa 



/l" 1 / 2 / dp 



-R r (-p)T r {- P y + T l {-p)R l (- P y 



+ h-^ 2 / dpT\-p)R l {-p)*i}(py qx/h e- lE v t/h 

J ~P0 

POO 

+ h~ 1 ' 2 d P T l {p)4,(p)e iqx/h e- lE » t/h , x > d. 



-iqx/h g—iE p t/h 



(21) 



Keep in mind that between — po < p < po, q is a positive imaginary number, q = i(\p 2 — 2mVo\) 1 ^ 2 . The second term 
in the third line and the fourth line come from a variable change p — > —p. 
Making use of Eqs. (fllHlft), there results the simple form 



(22) 



3 

/oo 
dpT l (p)^(p)e lqx/h e~ lE " t/R , x > d. 
-oo 

The equation can be put in an even more compact form as shown in the next section 
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IV. BOUND STATES 



None of the particles of the classical ensemble described in the introduction may be trapped by a potential well of 
the potential U(x). In quantum mechanics though, a bound state wave function extends exponentially beyond the 
potential limits, and may overlap with the initial state i^(0), even when this state is localized outside the potential 
limits. The contribution of these bound states to the wave function at t > is orthogonal to the scattering (continuum) 
part, and will remain spatially linked to the potential region at all times. 

A bound state with energy Ej corresponds to a simple pole of T l (p), or a zero of \/T l (p) on the positive imaginary 
axis, at pj = ijj, jj > 0. In this section we shall see that the bound state terms may be written as a residue 

<f> j (x)e- iE = t/h (E j \^(0)) = -2irih- 1/2 Res W l (p)^(p)e^ x/h e^ t/h ] , x >d (23) 

L J P = Pj 



where qj — i\ r~fj + Pq- The function ip(p) is defined on the complex plane by the integral 



^(P) = J^J2 j dx e- ipx/h ^(x, 0), (24) 

as has been used all along. 

For the formal treatment of bound states it is convenient to introduce Jost solutions /i (p, a:) and f 2 (p, x) of the 
Schrodinger equation. They are defined by the boundary conditions 

fi(p, x) — > e lqx / h , a,sx>d, and 

f 2 {p,x) - e -W*, x<c . (25) 

We generalize here the treatment of ref. Q for Vq — to the case Vq > 0. It is easy from (||) and ([?]) to obtain the 
explicit expressions of fi and f 2 at x < c and x > d respectively, 



iqx/h x> d 



fi(p,x) - I lpx/h ' j r - ivx /hR[(v) . „ ^ _ ( 26 ) 

e T --(p) +e T-( P ) ' a;>a (27) 
e -W» . x < c 

In particular, if for some value pj, 1/T l (pj) = 0, then f 2 and /i become real, proportional to each other, 

h = Ch, (28) 
and decay exponentially for x < c and x > d. Defining the normalization constant N by 



N 2 = / dx fl , (29) 



we may write the wave- function of the bound state \Ej) as 

4>i (x) = fi {p } ■ , x) = — f 2 (Pj,x). (30) 

The overlap between the bound state \Ej) and the initial state, which has support only on the interval [a, b] < c, is 
as usual 

(E^(0)) = f dx^(x)^(x,0). (31) 

On substituting the second equality of Eq. |3(]) in Eq. (|3l|), and taking into account definition ( p4| ) and Eq. (p7|), we 
obtain 

(E 3 \m) = ^ f dxf 2 ( Pj ,x)iP(x,0) = ^ f dxe-^/ n i;(x,0) = ^(Pj) ■ (32) 
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Therefore, for x > d, we have 



^(x)(E 3 \m) = (ihto,*)) mm) = (^#(o)) 



CN 2 



so that (^3|) is obtained. In the last line of ( |33| ) we have used || 



ResT l (p) p= 



p=Pi 



! -1 



P=Pi 



CN 2 



(33) 



(34) 



a relation that may be obtained by differentiating the Schrodinger equation for /i with respect to p, integrating /1/2 
between —R and R as R — > 00, and comparing with the derivative of the Wronskian of /1 and /2 with respect to p 
for p — pj, see also || . Finally, combining ( p2| ) and ( p3| ) we can write the transmitted wave packet very compactly as 

(35) 



i/>(x,t) = ^ 1/2 / dpT l (p)^(p)e lqx/h e- tE " t/n , x > d, 

JQ. 

where the contour Q goes from — 00 to 00 passing above the bound state poles, as shown in Fig. pi 




V. INITIAL STATE IN THE UPPER LEVEL 



Imagine now that the initial state has support only to the right, where the potential is Vb; that is, assume that the 
initial support of the state is the interval [a, b], such that d < a < b. In this case, instead of Eq. (19), we have 

T l (p)*$(q), p>p , 

<P + hK* = o)) - 1 T l ( P )y{-q), o<p< POl (36) 



(f) 1/2 {^(q)+R r (-p)*i>(-q) 



P < -Po , 



> 1/2 

where, as before, q — (p 2 — Pq) , with the branch cut going slightly below Im(p) = 0, and ip(q) is defined by Eq. 

©• 

Expression (|l^), as such, is also valid for x < c. Let us then substitute Eq. (|3^) in Eq. (|l8l), and use Eqs. (||) and 
(^). The result analogous to Eq. (|2l]) reads then, for x < c when the initial state is restricted at x > d, 



rl>(x,t) =J2 ( l>3(x)(E j m = 0))e- i 



Ejt/h 



G 



-I'd 



+ h~ 1/2 / dpT l {p)*^(q)e lpx/h e^ E " t/h 

J — OG 

+ h- 1/2 / dpT l {p)*^{q)e lpx/h e^ E " t/h 
Jpo 

rPo 

+ h~ 1 / 2 / d P T l { V )*^{-q)e lpx l h e- lE " t / h 



^R r ( P rT r (p)+T l ( P rR l (p) 



+ h- 1 ' 2 / dp 

J Vn 



Po 
po 



+ h- 1 ' 2 / dpT l {p)*R l (p)$(~q)e- ipx/R e- lE ? t/h , 
Jo 



-ipx /h —iEpt /h 



X < c. 



Making use of the unitarity relations, this expression can be written as 

^( s ,t) = J2<f>o( x )( E o\^(t = 0))e~ lE > t/h + h' 1/2 [ dpT l (~p)^(q)e lpx/h e- lE " t/n 



(37) 



(38) 



where the path of integration T goes from — oo to +00 immediately below the branch cut. Notice that for p = —17, 
with real positive 7, q becomes —i\f^ 2 + p\, which accounts for the minus sign in front of q in the third and fifth 
integrals of Eq. @. 

In order to include the bound states in a more compact expression, let us substitute the first equality of Eq. ( |30| ) 
in Eq. (EHI). Thus, using definition and Eq. (p6|), we obtain 



(39) 



(25^(0)) = i f dxf l (j> j ,x)Tp(x,0) = ± Ue^( I)0 ) = ^(- ft ), 



where, as before, qj = i\/J 2 + Pq- Therefore, for x < c, we have in this situation, using Eqs. (27) and (| 



^■(x)(^#(0)) = l—hfax)) (^|V(0)) = . 



= [ _ -ipjx/n 



(E 3 \m) 



h i ' 2 

CN 2 ' 



The pole of T'(p) at pj — ijj becomes a pole of T l (—p) at — pj. Thus, 



T\-p) = 



~ih 1 



P=-Pi 



CN 2 p + Pj ' 

and it follows that 

4> j (x)(E j \iP{t = 0)}e- lE J t/h = 2mh- 1/2 Res \T l (-p)e lpx/h iP(q)e-' lE ^ /h 

from which we can conclude that when the initial state is restricted to the interval [a, b], with d < a < b, 

if>{x, t) = h- 1 / 2 [ dpT l {-p)i}{q)e lpx/h e- lE ^/ n , 
iff 

for x < c, where f2' is a path of integration that goes from —00 to +00 below the branch cut and the poles —pj. 



(40) 



(41) 



(42) 



(43) 



VI. DISCUSSION 



The main results of this paper are Eqns. ( |35[ ) and (}43|), that provide simple, compact, and exact expressions for 
the transmitted wave packet in terms of the initial momentum distribution and bound state contributions. Whereas 
these expressions have been obtained for cut-off potentials, they may be easily generalized for potentials that decay 
fast enough, so that the initial wave packet does not overlap significantly with the potential region, and for x values 
where the asymptotic expressions for the scattering states are accurate. 
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These two equations make clear the need to include different contributions for the transmitted wave-packet. Even 
though tunnelling terms are of course included in theoretical analysis of wave packet collisions, the contribution of 
negative momenta is frequently overlooked in integral expressions of the wave function. Similarly, the contribution 
of the evanescent momentum region when the wave packet is initially on the potential upper level has been also 
disregarded. We have applied these equations in [|l0| to establish the relation between source boundary conditions, 
in which the wave function is specified at a point and for all times, with standard initial-value-problem boundary 
conditions where the wave function is specified in all space at t — 0. One further intended application is the study of 
arrival-time measurement models where a clock dial is coupled to the particle's motion in such a way that the particle's 
crossing stops the dial's motion Jll],[l2| • These models are described by step potentials of the form considered in section 
V. 

Let us now illustrate the formal results obtained above with some simple explicit computations. In all of them, the 
initial wave packet is taken as the ground state of an infinite well; it is located between a and 6, 

/ 2 \ 1/2 

ip{x, 0) = f — - J sin[(x-a)k w ]H(a,b), (44) 

where k w = it / {b — a), and 

H(a,b) = \l a<x<b, (45) 

v ' 10 otherwise. v ; 

See [|l0| for details of the analytical time dependence of ip(x,t). Here we shall evaluate the norm for x > 0, 

roc 

P T {t) = / dx\tp(x,t)\ 2 , (46) 



as well as the contributions from positive and negative momenta, the interference terms, and, whenever required, the 
contribution of bound states and evanescent waves. 
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FIG. 3. The contributions of positive and negative momentum to the probability of finding the particle in the half line x > 
for time t in the free case are shown with a dotted and a dashed line, respectively. The interference term is depicted with 
a dash-dotted line, and the total probability of finding the particle in the half line x > for time t is portrayed as a thick 
continuous line. The initial state is given by expression (fji|), with a = —2.01, b = —0.01. The mass m equals 1. Atomic units 
are used throughout. 



For instance, in the case of free motion, by defining the projectors 

/•OO 

P= dp\p)(p\ and 
Jo 

Q = i-P, 

Pt is decomposed into three terms, Pt = Pt,+ + Pt- + Pr.int- 



(47) 
(48) 
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Pr,-(t) = 

P T ,int{t) = 



dx\(x\Pil>(t))f 

> 

dx\{x\Q^{t))f 



dx2Rje[(x\Pip{t)}(Qip(t)\x)] . 



(49) 
(50) 
(51) 



These quantities are shown as functions of t for free motion in Fig. |^. The interference term is important only for 
early times, and the limit Pt(oo) = Pt,+ (oo) is clearly checked visually, although the decay of Pr,- is rather slow. 

In Fig. [| we show how the addition of an attractive delta potential, that creates a bound state, modifies the 
analogous contributions, since the bound state component has to be detracted from them. Pt, bound is constant in 
time, since there is only one bound state, and it is no longer true that the total final probability equals Pt,+ (oo); 



actually, we have Pr{oo) = P- a 



T, bound 



Pt. + {oo). In the picture we have only shown the interference term between 



positive and negative momenta, which is the most relevant one in this case. 



ra 




t (a.u.) 

FIG. 4. The initial state is as in Eq. (^), but multiplied by exp(ipa;), where p — 0.5 (all quantities in atomic units). Again 
we have a — —2.01 and b = —0.01. The motion of the particle is free everywhere, except for an attractive delta potential, 
U(x) = —UqS(x), with Uo = 1/8. The contribution of the bound state to the total probability of finding the particle in the 
half line x > is shown with a thin continuous line; the contributions of positive and negative momenta with a dotted and a 
dashed line, respectively. The interference term between positive and negative momenta is depicted with a dash-dotted line, 
and the total probability as a thick continuous line. 

In order to show the applicability of expression (|35|), we portray in Fig. || and Fig. ^| the squared moduli of the 
total amplitude and of different components of the amplitude for a step potential with an attractive delta potential; 
that is, the potential is V(x) = —Uq5(x) + Vq0(x). In this situation, additional to the positive and negative momenta 
components, we have evanescent waves and the contribution of the bound state. At the time of both pictures, (t = 5 
in atomic units), the evanescent waves are the major contributor to the wavefunction in the vicinity of the origin. It 
should be noticed though that interference terms are relevant in that region and elsewhere. 

The magnitude of the component of negative momenta, as well as that corresponding to the bound state, is rather 



small relative to the positive momenta component of the amplitude modulus squared, but it is nonetheless very 
important even for large distances because of the interference between components, as is clearly reflected in Fig. |5j. 
Fig. @ is added to show more clearly the relative importance of the negative momenta and bound state components. 
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FIG. 5. The particle of mass m = 1 moves in the potential V(x) — —UoS(x) + Vo9(x), where Uo = 1/8 and Vo = 1 (in 
atomic units). The initial state is as in Eq. (|l4|), with a = —2.01 and b = —0.01, but multiplied by exp(ipi), where p = 0.5. 
Time 4 = 5 has elapsed since the initial state was released. The moduli squares of different contributions to the amplitude 
are depicted: continuous thick line, total amplitude modulus squared; dotted line (resp. dashed line), square modulus of the 
contribution to the amplitude of positive (resp. negative) momenta; dot-dashed line, evanescent waves; continuous thin line, 
bound state contribution. 
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